f 


Historic,  archived  document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices 


i 


HELICOPTER  LOGGING: 

A  MODEL  FOR 
LOCATING  LANDINGS 


Louis  T.  Egging 
and 

David  F.  Gibson 


USDA  Forest  Service  Research  Paper  INT-155,  1974 
INTERMOUNTAIN  FOREST  &  RANGE 
EXPERIMENT  STATION 
Ogden,  Utah  84401 


USDA  Forest  Service 
Research  Paper  INT -155 
September  1974 


HELICOPTER  LOGGING: 
A  MODEL 
FOR  LOCATING  LANDINGS 


Louis  T.  Egging  and  David  F.  Gibson 


INTERMOUNTAIN  FOREST  AND  RANGE  EXPERIMENT  STATION 

Forest  Service 
U.  S.  Department  of  Agriculture 
Ogden,  Utah  84401 
Roger  R.  Bay,  Director 


THE  AUTHORS 


LOUIS  T.  EGGING  is  a  Graduate  Research  Assistant,  preseatl}' 
completing  work  for  his  master's  degree  at  Montana  State 
University,  Bozeman,  in  Industrial  and  Management  Engi- 
neering. He  received  hisB.A.  degree  in  Math  from  Carroll 
College,  Helena,  Montana,  and  his  B.S.  degree  in  Industrial 
and  Management  Engineering  from  M.  S.  U.  For  the  past 
3  years,  he  has  been  workit^  with  Intermountain  Station  in 
Bozeman  on  Logging  Systems  studies. 

DAVID  F.  GIBSON  is  Associate  Professor,  Industrial  and  Man- 
agement Engineering,  Montana  State  University.  For  the 
past  2  years,  he  has  been  working  cooperatively  with  the 
Intermountain  Station's  Forest  Engineering  Research  Unit 
at  Bozeman  as  a  Research  Industrial  Engineer  under  the 
Intergovernmental  Personnel  Act.  He  holds  B.S.,  M.S., 
and  Ph.  D.  degrees  in  Industrial  Engineering  from  Purdue 
University. 


CONTENTS 


Page 


INTRODUCTION   1 

STATEMENT  OF  THE  PROBLEM   2 

COMPUTER  MODEL   4 

Solution  Procedure   4 

Program  Input   7 

Examples   11 

Problem  1   11 

Problem  2   19 

Limitations   26 


LITERATURE  CITED 


27 


ABSTRACT 


Presented  are  a  model  and  an  accompanying  computer  program 
that  optimally  locate  landing  areas  for  a  helicopter  logging  operation. 
Given  a  haul  road,  unit  centroids,  volumes  of  timber  to  be  harvested, 
and  helicopter  operating  parameters,  landings  are  located  so  as  to 
minimize  yardit^,  hauling,  and  landing  construction  costs.  The 
model  considers  constraints  such  as  areas  that  are  not  suitable  for 
landings  and  topographical  obstacles.  Written  in  FORTRAN  IV,  the 
computer  program  affords  several  evaluation  and  output  options. 
Two  examples  are  provided. 


INTRODUCTION 


In  recent  years,  helicopter  logging  has  emerged  from  experimental  stages  to 
become  a  viable  alternative  to  conventional  harvesting  systems.    Helicopter  logging 
has  a  number  of  attractive  advantages  under  certain  conditions  and  management  objec- 
tives.   These  include  minimal  environmental  impact  and  the  ability  to  harvest  hereto- 
fore inaccessible  timber.     As  a  result,  the  volume  of  timber  harvested  by  means  of 
helicopter  systems  has  been  increasing  steadily  each  year. 

^     Binkley  [1972]  Edholm  (1973,  1974),  and  Stevens  (1972,  1973)  have  reported  the 
general  production  and  cost  aspects  of  helicopter  logging.    The  studies  show  that  heli- 
copter logging  is  expensive  and  warrants  efforts  to  design  an  operating  system  in  the 
most  economical  manner  possible.     Gibson  (1974)  developed  a  cost  model  that,  given 
parameters  of  a  system,  specified  optimum  refueling  of  the  aircraft. 

This  paper  presents  another  tool  for  designing  and  managing  a  helicopter  logging 
sale  by  presenting  a  computer  model  that  specifies  the  number  and  location  of  landings 
to  which  the  timber  should  be  yarded.     Considerations  include  areas  not  suitable  for 
landings,  and  flight  path  restrictions  such  as  ridges  or  other  obstacles.  Alternative 
landing  sites  are  evaluated  by  total  cost  of  yarding,  landing  construction,  and  hauling. 


STATEMENT  OF  THE  PROBLEM 


Timber  felled  and  "marked  for  turns"  (logs  identified  to  be  grouped  together  as  a 
load)  is  yarded  to  landings  that  are  located  on  or  in  close  proximity  to  a  preexisting 
road.     The  timber  is  then  loaded  on  trucks  and  hauled  to  the  mill.     Location  and  the 
number  of  landings  affect  (1)  yarding  costs,   (2)   landing  construction  costs,  and 
(3)  hauling  costs.    Topographical  features  often  influence  flight  paths  and  landing 
selection.     This  paper  is  addressed  to  the  problem  of  locating  landings  to  minimize 
costs  and  yet  satisfy  constraints. 

Locating  helicopter  landings  is  a  special  case  of  the  generalized  Weber  problem, 
which  is  concerned  with  locating  facilities  within  a  given  space  so  as  to  minimize  a 
cost  function.     Many  variations  of  the  problem  have  been  defined  and  include  the  con- 
sideration of  various  measures  of  distance  (rectilinear,  Euclidean,  square  of 
Euclidean),  costs,  and  other  features.     Notable  treatments  of  the  problem  include  those 
by  Kuhn  and  Kuenne  (1962),  Bellman  (1965),  Cooper  (1963,  1968),  Hakimi  (1964),  Levy 
(1967),  and  Cabot  and  others  (1970). 

The  treatment  presented  in  this  paper  is  unique  in  several  respects.  First, 
although  Euclidean  distances  are  used  to  determine  yarding  costs,  obstacles  may  be 
present:  a  flight  path  may  be  a  series  of  linear  line  segments  routed  around  a  ridge 
or  other  obstacle,  as  opposed  to  a  single  path.    Secondly,  the  costs  of  constructing 
landings  has  not  only  a  fixed  component,  but  also  a  variable  one  as  well,  which  is 
dependent  upon  the  amount  of  timber  yarded  to  it.     Also,  a  second  transportation  cost-- 
the  hauling  cost--is  included. 

The  extensive  literature  related  to  the  generalized  Weber  problem  will  not  be  re- 
viewed here.     Nor  will  the  computational  techniques  previously  developed  be  enriched 
to  extend  the  model  to  more  general  applications.     This  may  be  done  in  a  subsequent 
paper.     The  purpose  of  this  paper  is  to  define  the  problem  of  designing  and  managing 
helicopter  logging  systems  and  to  provide  a  FORTRAN  IV  program  for  locating  landings. 
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A  statement  of  the  problem  follows: 


Given : 


1.  a  haul  road  that  can  be  represented  by  a  set  of  linear  line  segments, 

2.  units  to  be  logged,  together  with  their  respective  centroids  and  timber 
volumes , 

3.  road  segments  infeasible  for  landing  locations, 

4.  ridges  or  other  topographical  obstacles  (which  can  be  represented  by  linear 
line  segments)  that  may  preclude  direct  flight  paths  between  units  and  land- 
ings, and 

5.  system  costs  and  operating  characteristics. 


find: 


1.  the  number  of  landings, 

2.  the  location  of  landings  on  the  haul  road,  and 

3.  the  allocation  of  logging  units  to  landings, 

so  as  to  minimize  the  sum  of: 

1.  yarding  cost, 

2.  landing  construction  cost,  and 

3.  hauling  cost. 

Figure  1  illustrates  a  helicopter  logging  layout  with  four  units  and  a  haul  road 
represented  by  eight  linear  line  segments.    On  two  of  the  units,  ridges  may  prevent 
direct  flights  to  landings.     The  problem  is  to  find  the  number  (1,  2,  3,  or  4)  and 
location  of  landings  and  the  allocation  of  logging  units  to  landings  that  will  minimize 
costs . 


Figure  1. — A  heliooptei' 
logging  layout. 
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COMPUTER  MODEL 


Solution  Procedure 

A  simplified  flow  chart  of  the  landing  location  model  is  given  in  figure  2.  Block 
1  represents  the  input  required  for  solution  of  the  problem.     This  includes  information 
concerning  the  units  to  be  logged,  the  haul  road,  operating  characteristics  of  the 
logging  system,  landing  construction,  hauling  and  operating  costs,  and  topographic 
features  of  the  area. 

Block  2  represents  the  establishment  of  a  scheme  for  allocating  logging  units  to 
landings.     For  example,  if  three  units  are  to  be  logged,  timber  in  the  first  unit  may 
be  yarded  to  one  landing  and  timber  in  the  other  two  units  may  be  yarded  to  a  second 
landing.    Another  allocation  scheme  would  be  to  yard  all  timber  to  one  landing.  The 
program  has  the  facility  to  search  all  possible  allocation  schemes  or  only  a  particular 
subset  that  the  user  is  especially  interested  in. 

Once  an  allocation  scheme  is  fixed,  the  program  then  finds  the  best  location  for 
the  specified  landings.    The  program  can  search  along  the  entire  road  for  these  loca- 
tions or  will  search  only  specified  points.     In  either  case,  it  begins  the  search  by 
fixing  a  location.     This  operation  is  represented  by  block  3. 

Next,   flight  paths  and  distances  between  the  units  and  the  landings  are  deter- 
mined.   When  direct  flight  paths  are  precluded  because  of  ridges  or  other  obstacles, 
the  flight  path  is  determined  to  be  the  shortest  distance  composed  of  two  linear  line 
segments  around  the  obstacle,  as  shown  in  figure  3.     (This  operation  is  represented  by 
block  4  in  figure  2.).     Euclidean  distances   (including  elevation)  are  used.     In  a 
situation  similar  to  that  shown  in  figure  3,  since  the  centroid  is  used  in  describing 
the  average  flight  path,  all  timber  in  the  xmit  is  routed  around  the  ridge.     If  the 
ridge  runs  through  a  unit  between  the  road  and  the  centroid,  the  analyst  may  wish  to 
define  two  units,  one  on  each  side  of  the  ridge. 
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FIX  ALLOCATION 
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FIX  LANDING 
LOCATION  3 


DETERMINE 
FLIGHT  PATHS  4 


CALCULATE 
TOTAL  COST  5 


SELECT  OPTIMUM 
SOLUTION  - 

y 

\ 

OUTPUT 

\  w 

Figure  2. — Simplified  flew  chart  of  landing  location  model. 
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UNIT  CENTROID 


LANDING 


Figicce  S. — Helioopter 
flight  path  around 
ohstaale. 


Total  costs  for  the  allocation  scheme  and  location  under  consideration  are  then 
calculated  as  represented  by  block  5.     Included  are  yarding  costs,  landing  construction 
costs,  and  hauling  costs.    Abbreviated  formulas  for  each  of  these  are  given  in  figure  4. 
Yarding  costs  include  the  consideration  of  the  average  speed  of  the  helicopter  when 
loaded  and  the  average  speed  when  empty.     Landing  construction  costs  can  vary  with 
location  and  have  both  fixed  and  variable  (with  weight  yarded)  components. 

Block  6  of  figure  2  signifies  the  program's  check  to  see  if  all  locations  have 
been  evaluated;  if  they  have  not,  the  cost  of  the  next  location  is  calculated.  The 
next  location  can  be  at  a  specified  increment  along  the  road  or  at  a  particular  loca- 
tion specified  by  the  user.    This  process  is  continued  until  all  locations  have  been 
evaluated  for  a  given  allocation  scheme. 


YARDING  COST  =     I     (DISTANCE  ^  SPEED  X  ^^^^X  #  TR  IPS) 
UNITS  '  "^"^ 


LANDING  CONSTRUCTION  COST  =        S        (FIXED  COST  +  VARIABLE  COST) 

LANDINGS 

HAULING  COST  =        S        (DISTANCE  X  #  TRIPS  X  nicxf lire) 
LANDINGS  UlblAlMUb 


Figure  4. — Abbreviated  oost  equations- 
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Intermediate  output  of  suboptimum  solutions  can  be  obtained  as  shown  by  block  7. 
That  is,  the  program  can,  at  the  discretion  of  the  user,  print  out  the  best  solution 
for  a  given  allocation  scheme. 

Once  all  locations  have  been  evaluated  for  the  current  allocation  and  the  minimum 
cost  location  has  been  selected,  the  program  checks  to  see  if  all  allocation  schemes 
have  been  examined  (block  8).     If  all  schemes  have  not  been  evaluated,  the  process 
cycles  back  to  block  2;  otherwise  the  optimal  solution  is  selected  from  previous  cal- 
culations (block  9).     Block  10  represents  final  output,  which  may  be  in  several  forms 
at  the  option  of  the  user,  including  graphical  plotting. 

Program  Input 

The  computer  program,  written  in  FORTRAN  IV,  is  available  on  request  from  the 
Forestry  Sciences  Laboratory  (ATTN:  David  F.  Gibson),  USDA  Forest  Service,  Inter- 
mountain  Forest  and  Range  Experiment  Station,  P.  0.  Box  1376,  Bozeman,  Montana  59715, 
Input  to  the  program  is  documented  via  comments;  however,  a  more  thorough  explanation 
follows : 


Card  Type      Card  Col.      Variable  Explanation 

Card  Set  1:  Parameter  Cards 
1 


Title  Card  1-48 


ITITLE      Any  48  alpha-numeric  characters  can  be  used  to  de- 
fine the  title  of  a  particular  run.    This  title 
will  be  printed  as  a  heading  on  output  from  the 
line  printer  and  as  a  caption  on  plots. 


Data  Card  1-2 
3-4 

5-6 

7-16 
17-26 
27-36 

37-46 

47-56 


NUNIT       The  number  of  units  that  are  to  be  logged.  (Maxi- 
mum of  8.) 

NROAD       The  haul  road  is  defined  as  a  set  of  linear  line  seg- 
ments.   NROAD  is  the  number  of  linear  line  segments 
representing  the  road.     (Maximum  of  20.) 

NRWL         Certain  segments  of  the  road  can  be  defined  as  in- 
feasible  locations  for  landings.     NRWL  is  the 
number  of  such  restricted  segments.     (Maximum  of  10.) 

P  Average  weight  per  load  of  the  helicopter. 

OC  System  operating  cost  expressed  in  $/min. 

SI  Average  speed  of  helicopter  when  loaded  expressed 

in  ft/min. 

SO  Average  speed  of  helicopter  when  unloaded  expressed 

in  ft/min. 

CPM  Cost  to  haul  logs  from  landing  to  the  mill  expressed 

in  $/mile. 


57-66 


WPL 


Average  weight  of  truckload  expressed  in  lb. 
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Card  Type      Card  Col.  Variable 


Explanation 


67 -lb  XINC  XINC  is  the  increment  the  program  advances  along 

the  road  in  search  of  the  optimum  locations. 
(Expressed  in  feet.) 

77  TOUT         A  variable  to  control  the  amount  of  output  of  the 

program.     When  TOUT  =  0,  the  complete  output  is 
printed,  which  includes  input  values,  the  solution 
for  each  allocation  scheme,  and  the  optimal  solu- 
tion.    If  TOUT  is  set  equal  to  1,  a  printout  of 
the  input  values  and  the  optimal  solution  is  made. 
A  printout  of  only  the  optimal  solution  is  given 
when  lOUT  =  2. 

78  IPLOT       This  variable  determines  whether  the  optimal  solu- 

tion is  drawn  graphically  by  the  CALCOMP  plotter. 
If  IPLOT  =  1,  a  plot  is  made.     If  no  plot  is  de- 
sired or  a  CALCOMP  plotter  is  not  available,  set 
IPLOT  =  0. 

79-80  NALT         Specific  landing  locations  can  be  investigated 

using  this  variable.     If  NALT  =  0,  the  road  is 
searched  every  XINC  ft  for  the  best  location. 
Otherwise  NALT  is  the  number  of  alternative  landing 
locations  to  be  tested. 

3  This  card  is  required  only  if  IPLOT  on  card  2  is  1. 

Plot  Card      1-2  NCOP         The  number  of  copies  of  each  plot  desired. 

3-7  PH  The  height  in  inches  of  the  plot. 

8-12  PW  The  width  in  inches  of  the  plot. 

Card  Set  2:  Unit  and  Ridge  Cards.     This  card  set  contains  unit  and  ridge  constraint 
cards.    There  is  one  card  for  each  unit.    Therefore,  there  are  I  =  1,  2,  NUNIT 
unit  cards  in  this  set.     In  addition,  there  is  one  card  for  each  ridge  constraint,  but 
a  maximum  of  one  ridge  constraint  per  unit.     Ridge  constraint  cards  follow  the  unit 
cards  for  the  unit  to  which  they  pertain. 


Unit  Card 


1- 

10 

UX(I) 

X  -  coordinate  of  unit  I. 

11- 

20 

UY(I) 

Y  -  coordinate  of  unit  I . 

21- 

30 

UZ(I) 

Z  -  coordinate  of  unit  I  (elevat 

ion)  . 

31- 

40 

wci) 

Total  weight  of  timber  in  unit  I 

expressed  in  lb 

41 

IRC 

If  unit  I  has  a  ridge  constraint 

associated  with 

IRC  is  input  as  an  R  and  then  this  card  is  immedi- 
ately followed  by  a  ridge  constraint  card. 
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Card  Type      Card  Col.      Variable  Explanation   

2  This  card  type  is  required  only  if  an  R  appears  in  column  41  of  the  pre- 

ceding card. 


Ridge  Con- 

1-10 

RCX(I ,1) 

X  - 

coordinate 

of 

beginning  of  ridge  constraint  on 

straint 

unit  I . 

Card 

11-20 

RCY(I,1) 

Y  - 

coordinate 

of 

beginning  of  ridge  constraint  on 

unit  I . 

21-30 

RCZ(I,1) 

Z  - 

coordinate 

of 

beginning  of  ridge  constraint  on 

unit  I . 

31-40 

RCX(I,2) 

X  - 

coordinate 

of 

end  of  ridge  constraint  on  unit  I 

41-50 

RCY(I,2) 

Y  - 

coordinate 

of 

end  of  ridge  constraint  on  unit  I 

51-60 

RCZ(I,2) 

Z  - 

coordinate 

of 

end  of  ridge  constraint  on  unit  I 

Card  Set  3: 

Roadway  Cards.  This 

set 

of  cards  contains  data  pertaining  to  the  linear 

line  segments  that  define  the  road.  This  set  is  not  required  if  NROAD  (card  set  1, 
card  type  2,  columns  3-4)  is  0.    That  is,  the  program  can  be  run  without  defining  a 


road.  In  such  a  case,  NALT  would  be  greater  than  0,  and  the  analyst  would  be  investi 
gating  a  particular  set  of  landing  location  alternatives. 

1 


Origin  Card 

1-10 

RX(0) 

X  -  coordinate  of  beginning  point  of  road. 

11-20 

RYCO) 

Y  -  coordinate  of  beginning  point  of  road. 

21-30 

RZ(0) 

Z  -  coordinate  or  i^oginning  point  of  road. 

2 

There  must 

be  one 

of  these  cards  for  each  linear  line  ^^^^^^^ 

Segment 

1-10 

RX(I) 

X  -  coordinate  of  ending  point  of  road  segment  i. 

Cards 

11-20 

RY(I) 

Y  -  coordinate  of  ending  point  of  road  segment  I. 

21-30 

RZ(I) 

Z  -  coordinate  of  ending  point  of  road  segment  I . 

31-40 

FC(I) 

Fixed  cost  of  constructing  a  landing  on  road  segment 

expressed  in  $. 

41-50  VC(I)        If  the  cost  of  a  landing  on  road  segment  I  is  also 

dependent  upon  the  amount  of  timber  yarded  to  the 
landing,  then  the  variable  VC(I)  can  be  employed  to 
express  this  cost  in  terms  of  $/lb. 
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Card  Type      Card  Col.      Variable  Explanation  

Card  Set  4:     Restricted  Segment  Cards.    This  set  of  cards  defines  those  segments  of 
the  road  where  it  is  not  possible  to  construct  a  landing.     IF  NRWL  (card  set  1,  card 
type  2,  columns  5-6)  is  0,  this  set  is  not  required.    Otherwise  there  will  be  NRWL 
cards  in  this  set,  one  for  each  restricted  segment. 

1 


Restricted 

1- 

10 

.  RXLCI,1) 

X  - 

coordinate 

of 

beginning  of  restricted  segment 

Segment 

Cards 

11- 

20 

RLY(I,1) 

Y  - 

coordinate 

of 

beginning  of  restricted  segment 

21- 

30 

RLX(I,2) 

X  - 

coordinate 

of 

end  of  restricted  segment. 

31- 

40 

RLY(1 ,2) 

Y  - 

coordinate 

of 

end  of  restricted  segment. 

Card  Set  5:    Alternative  Landing  Cards.    The  card  set  contains  information  about  spe- 
cific landing  location  alternatives  the  user  may  want  to  investigate.    This  set  is 
not  required  if  NALT  (card  set  1,  card  type  2,  columns  79-80)  is  0.     Otherwise  there 
will  be  NALT  cards  in  this  set,  one  for  each  alternative  landing. 

1 

X  -  coordinate  of  landing  I. 
Y  -  coordinate  of  landing  I. 
Z  -  coordinate  of  landing  I.- ' 

The  fixed  cost  of  constructine  a  landing  ar  location  I 
expre'^s'^''  ^"  ^' 

The  variable  (with  amount  of  timber  yarded  to  the 
landing)  cost  of  constructing  a  landing  at  loca:tion  1 
expressed  in  $/lb. 

Card  Set  6:     Allocation  Scheme  Cards.    This  card  set  may  be  repeated  as  many  times  as 
desired  by  the  user,  depending  upon  the  number  of  allocation  schemes  to  be  investi- 
gated.    A  new  set  is  required  for  each  different  number  of  landings  investigated. 

1 

Allocation      1  ID  The  variable  ID  can  be  set  equal  to  zero  (blank)  or 

Control  Card  input  as  A,     If  ID  =  0,  then  the  second  card  type 

will  contain  allocation  matrix  cards  and  the  program 
.    -      will  search  all  possible  locations  on  the  road  for 
each  allocation  combination  and  minimize  the  total 
cost  in  deriving  the  best  solution.     If  ID  =  A,  the 
second  card  type  will  contain  alternative  allocation 
cards,  and  the  program  will  compare  and  print  the 
cost  of  specific  location-allocation  schemes. 


Alternative    1-10  ALTX(I) 
Landing 

Card  11-20  ALTY(I) 

21-30  ALTZ(I) 

31-40  FC(I) 


2-3  NLAND       The  number  of  landings  to  be  employed  in  the  alloca- 

tion scheme (s)  defined  in  the  second  card  type. 

4-5  NCOM         The  number  of  allocation  schemes  to  follow  in  the 

second  card  type. 
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Card  Type      Card  Col.  Variable 


Explanation 


This  type  of  card  is  used  only  if  ID  =  0  (or  blank) . 
If  ID  =  0  (or  blank),  there  will  be  NCOM  cards  of 
this  type.    Each  allocation  scheme  will  have  a  card 
whose  entries  are  the  elements  of  an  allocation 
matrix.     Rows  represent  units  and  columns  represent 
landings.     If  unit  I  is  to  be  assigned  to  landing  J, 
then  the  I, J  element  of  the  matrix  is  1,  otherwise 
it  is  0.    Entries  are  read  in  by  column  with  each 
entry  allotted  one  card  column.    These  entries  can 
perhaps  be  best  understood  by  considering  the  examples 
presented  later  in  the  text  and  the  description  given 
via  comment  cards  in  the  program. 


This  type  of  card  is  used  only  if  ID  =  A.     If  ID  =  A, 
there  will  be  NCOM  cards  of  this  type.    Each  allo- 
cation combination  will  have  a  card  which  will  have 
NUNIT  (number  of  units  to  be  harvested)  entries.  The 
entry  in  the  Ith  column  will  be  the  alternative  land- 
ing number  to  which  unit  I  is  to  be  yarded.  These 
entries  can  be  perhaps  best  understood  by  considering 
the  examples  presented. 

Card  Set  7:     End  of  Job.    The  last  card  set  contains  one  card,  a  blank  card,  which 
signifies  that  the  problem  input  is  complete.     Problems  can  be  run  "back  to  back"  by 
repeating  card  sets  1  through  7  successively. 


PToblem  1 

As  an  example  of  the  use  of  the  model,  consider  the  problem  portrayed  in  figure  1. 
Figure  5  gives  the  layout  with  coordinates,  elevations,  and  other  data  required  for 
solution.     Input  to  the  problem  as  prescribed  in  the  preceding  section  is  illustrated 
on  a  FORTRAN  coding  sheet  in  figure  6.    Output  from  the  program  on  the  line  printer 
is  given  as  figure  7.     Since  the  input  variable  "lOUT"  (column  77,  card  2)  was  set 
equal  to  "0,"  complete  output  was  obtained.    That  is,  input  data,  intermediate  solu- 
tions, and  the  optimal  solution  were  printed  out.    The  solution  for  case  7,  when  two 
landings  are  considered,  indicates  that  both  landings  are  in  the  same  location. 
Since  two  landings  were  specified  in  the  allocation  scheme,  the  fixed  cost  is  added 
twice.    However,  this  solution,  relative  to  location,  should  then  be  the  same  as  for 
the  allocation  scheme  with  one  landing.     Indeed  this  is  the  case  as  indicated  in  the 
output.    The  output  obtained  from  the  plotter  (obtained  because  input  variable  IPLOT, 
column  78,  card  2,  was  set  equal  to  1)  is  shown  in  figure  8. 


Al location 
Matric 
Cards 


Alternative 
Allocation 
Cards 
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8,000  ■■ 


(6,900;  8,300;  800) 
/ 


7,000-- 


6,000  ■  ■ 


5,000-- 


4,000  ■■ 


3,000- 


2,000  ■  - 


UN  IT  4  (5,500;  7,500;  900)  , 

(3,800;  7,300;  1,000)  •  / 

1,700,000 


'  (6,000;  7,400;  700) 


(4.500;  6,000;  800) 


UNIT  3 

(1,600,  4,200;  800)  COORDINATES 
950,000  UNIT  WEIGHT 


(3,000;  3,500;  700) 
^\OGt  ^  .  -  '  '  * 


(5,600;  5,700;  600) 


(5,300;  4,200;  500) 


(1,000;  3,000/ eoo; 


(4,500;  3,600;  400) 
(3,900;  3,100;  300) 


UNIT  2 
(6,800;  3,200;  600) 
1,200,000 


RESTRICTED  ROADWAY 


X     (1,400;  2,000;  200) 


1  OOOf        /  <500;  1,100;  100)  UNIT  1 

/  (2,500;  1,200;  500) 

/  800,000 


1,000 


2,000  3,000 


4.000 


5,000  6,000 


7,000 


Figure  5. — Layout  of  helicopter  logging  problem  1  (from  figure  1) . 
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IBM 


FORTHAN  Coding  Fortn 


GX28-7327-6U/M050' 
Printed  in  U  S  A 


I 

FORTRAN  STATEMENT 

E 

xaMp 

_ 

E 

PR 

©BLEW 

»»'1 

[ 

I'll 

4 

8 

2 

8000 

■  '  'V 

8 

[3  0 

olo 

if 

4'obo 

.4 

0  ■ 

60000 

1  UU  • U  T  U 

r  i 

1 

1 

.  :  18 

•■5  ■ 

2 

500  . 

1  200  . 

500  . 

80 

0000  . 

6 

800  . 

3200  . 

600  . 

1  20 

0000  . 

1 

600  . 

4200  . 

800  . 

95 

0000  . 

R 

-J  

1 

ooo . 

3000!. 

&0l0 

1 

300  a. 

3  50'0. 

7  0:0. 

3 

8iOO. 

7  30  0  . 

1  OiOiO 

1 

70 

000  0 . 

R 

4 

5.0;O  . 

■ 

60O'0  . 

800 

5510  O: . 

7iS00;. 

,9i00Lj 

0  • 

0  . 

0  . 

500  . 

1  1  00  . 

100  . 

1  500  . 

.  0001 

1 

400  . 

2000  . 

200  . 

1  800  . 

.  0001 

1 

3 

gio'oi . 

31!0!0i. 

3'OlOi . 

1  60  Oj . 

.  0:0  0  2 

4 

soo 

3  600'  • 

400 

1 ,40  01  • 

.0 

I 

5 

3:00 

4  200. 

5.0O 

1  lOO!. 

.  00:0  3 

5 

600  . 

5700  . 

600  . 

1  300  . 

.  0002 

6 

000  . 

7400  . 

700  . 

1  300  . 

.  0003 

5 

900  . 

8300  . 

800  . 

1  300  . 

.  0003 

4- 

1 

4O0: . 

2  OOOl . 

3'9'0l01 . 

3  1O  0' . 

i 

1 

4- 

1 — 

4 

5I0;0| . 

3:6:O0  . 

i 

SGO'O  ■ 

■ 

14 

1 

— 

i  : 

i 

1 

0000 

1 

0000 

1  0000 

1 

3  6 

1 

1000 

0 

1  000 

01 

IBM 


FORTRAN  Coding  Form 


I 

FORTRAN  STATEMENT 

1 

0 1 :00 

1 

0,000 

01 

] 

1 

-4- 

-idU- 

1 

- 

i  1 

1  :  '  : 

1 

00:1  0 

1 

0  00  0 

1:0 

1 

1 

4- 

:  1  1 

■  '  1  — 

0 

1  1:01 

0 

0000 

01 

— M 

I  1  j 

0 

1011 

0 

0000 

1  0 

0 

0111 

0 

000  1 

00 

2  7 

1 

l'l  00 

0 

01 

1 

0110 

1 

00 

0 

11  1  1 

0 

00 

1 

1010 

0 

1  0 

1 

1000 

0 

1  1 

1 

01  00 

1 

01 

1 

001  0 

1 

1  0 

1  1 

1 

1  11 

0 

000 

Figure  6. — Input  data^  helicopter  logging  problem  1. 
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Figure  7.  — Line  printev  output. 


helicopter  togging  problem  1. 


MKlT  NRC*C   t.«WL       P  OC  SI  so  CPM 

*  8  2  •ecOCe  Ot      .IgOOE  O?      •SOOOF  .^OOCE  Ot      .itOOCE  00 

V.PI  XiNf  lOlT  IPLTT  naLT    KCCP     Ph  PV, 

.600CE   C5      •ICOOF  03     C         1  0  1      lltOC  8.50 


ttft   LMT   AND   RIDGE   CONSTRAINT  CARDS 


I 

I  y  ( I ) 

L  Y  (  I  ) 

lj'Z(  I  ) 

W  (  I  ) 

IRC 

RCX(  I, 

1  ( 

1 

.?5CF 

Ck 

.i2cr 

04 

.500E 

03 

«  flCOE 

06 

.  COCE 

00 

«  f  80f 

Ct 

.3?0E 

04 

.  600E 

03 

•  1?0F 

07 

•  COOE 

on 

3 

04 

.-♦gOE 

04 

.800F 

03 

•9!:0F 

06 

R 

•  ICCE 

04 

4 

•  3PCF 

Ct 

.73CE 

04 

•  lOOE 

04 

•  170F 

07 

R 

.45CE 

04 

RCY (  III) 
•COCE  00 
•COCE  00 
.30CF  04 
•fOCE  04 


RC7( I<1 ) 
«00CE  CO 
•COCE  00 
.600E  03 
.80CE  03 


RCX( 1,2) 
•COOE  00 
•COCE  00 
•3C0E  04 
•550E  04 


RCY( I<2) 
•COOE  OC 
•OOCE  OC 
.35CE  04 
.75CE  04 


RC7(1,2) 
•COOE  00 
•COOE  00 
•70CE  03 
•?OCE  05 


•  RCaCuAY  CARDS 


I 

RX  (  1  ) 

RY(  I  ) 

RZ( I ) i 

rc(  I ) 

vC  (  I  ) 

c 

•  CCOF 

CC 

•  COOE 

OC 

•  OOOE 

00 

1 

•  5C0F 

Or? 

•  llOE 

04 

•  lOOE 

03 

•  150F 

04 

.1 COE- 

03 

2 

•  1  40F 

C4 

.2C0E 

04 

.200E 

03 

•  1  ROE 

04 

.1C0E- 

03 

3 

•  390F 

04 

.310E 

C4 

.300E 

03 

•  IfeOE 

04 

•PCCE- 

03 

4 

•  450F 

04 

.360E 

C4 

.400E 

03 

•  140F 

04 

•  OOCE 

CO 

5 

•  E3CE 

C4 

.420E 

04 

•  SOOE 

03 

•  110E 

04 

•  3CCE- 

03 

6 

•  560F 

04 

.570E 

04 

•  600E 

03 

•  130F 

C4 

.2C0E- 

03 

7 

•  600F 

C4 

.74CE 

04 

.700E 

03 

•  130r 

C4 

.300E- 

03 

8 

•  690F 

C4 

•  830E 

04 

•  SOOE 

03 

•  130E 

04 

.3C0E- 

03 

RFSTRICTED 

RC ADwAY 

SEGMENTS  ^ 

»»»» 

I 

RLX (  1/ 

1  ) 

RLY (  I  i 

1  ) 

RLXt  I, 

2) 

RLY(I>2) 

1 

•  1  40F 

C4 

.20GE 

04 

•  390E 

•  310F 

04 

2 

•  450E 

C4 

.360E 

04 

•  560E 

0* 

•  570F 

04 

NU^'8ER  OF  LMTS  -  4 
^U^'BER     OF   LANDINGS  ■  4 


c*sE  #  i: 


LANDING     C00R0lNATES(X<yi7)        UNIT  DISTANCE 


#1  (  1400*  2000,  200) 

#2  I  4500'  3600,  400) 

#3  (  3900,  3100,  300) 

#4  (  5989'  7354'  697  1 


YARDING  COST 


13443^47 


HAULING  COST   ■  « 


39.89 


TCTaL  CCST 


19913^36 


f  1 
#2 
#3 
#4 


1392  FT 

2343  FT 

2598  FT 

2263  FT 


LANCING  BUILDING  COsT  ■  «  6429.99 
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NUf'BPR   OF   LMTS   -  4 
NUfBER     CF   LANDINGS  -  3 


LANCING 

COORDINATE'S  (  X#  Yi  7  ) 

LN  IT 

C ISTANC 

#1  ( 

hi^37'  354fl/ 

389  t 

#1 

3C46  F 

#2 

2397  F 

#?  ( 

3900<  3100< 

30r ) 

#3 

2598  F 

«3  1 

5989*  7354/ 

697  ) 

«% 

2263  F 

YaRC ING 

COST   -   •  \^Zek.70 

Lancing 

building  cost 

•  • 

461C 

.CC 

Hall ING 

COST   -  • 

43.?'t 

TCTaL   CCST   -  $     19917. 94 


r  A  s  F  #    2  : 


LANDING 

C00RDlNATES(X*Vi7l 

UNIT 

DISTANC 

i»l  ( 

1400*  2000* 

200  ) 

*1 

1392  F 

#3 

2505  F 

»2  ( 

4500*  3600* 

400  ) 

#2 

2343  F 

#3  ( 

5989*  7354* 

697  ) 

1114 

2263  F 

YARDING 

COST  ■  •  13327.89 

Landing 

BUILDING  COST 

■  c 

5220 

.CO 

HALL ING 

COST   -  $ 

36  .61 

TCTaL   CCST   ■  $  1R58'».50 


LANDING     COORDINATFS ( X* Y, 7  1 


#1 


(    4437*    3548*  389) 


«2  (    4500*    3600*  4001 

»3  1    3900'    3100*  3001 

YARDING   COST   ■   t  18728.51 
Lancing  BUILDING  COST  ■  • 
HALl  ING   COST   ■   •  34 .08 

TCTaL   CCST  -  «  23022.59 


UNIT  DISTANCE 


#1 
«4 
#2 
«3 


3046  FT 

3854  FT 

2343  FT 

2S98  FT 


4260. CC 


CASE   *  4: 


LANDING 
*1 


C0CRDINATFS(X*Y.7I 
(    4437*    3548*  3891 


#2  (    1400*    2000*  200) 

*3  (    5989*    7354*  697) 

YARDING  COST  -  «  13955.30 
Lancing  BUILDING  COST  >  • 
HALLING   COST   •   «  40.61 
TCTaL   cost   ■  C  18965.91 


LMT  DISTANCE 
♦2       2397  FT 
2940  FT 
1392  FT 
2263  FT 


«3 
#1 
*4 


4970. CC 


CASE    »      5  1 


LANDING     COORDINATES*  X»Yi7  )  LMT  DISTANCE 

*1          (    5600*   5700*      600)  #2       2773  FT 

#4       2647  FT 

»2          (    1400*    2000*      200)  »1       1392  FT 

#3  (  3900*  3100*  300)  #3  2598  FT 
YARDING   COST   •   •  14976.89 

Landing  buildinq  cost  ■  «  5330. CO 

HALLING   COST   ■   t  40.08 
TCTaL   COST   -  •  20346.97 


15 


CASE  »  6: 


LANniNG     CCORDP.ATFS(      Yj7  I 
#1  (    560C*    570C»  60C) 

#?       (  140C'  ?ooo»  ?or) 

#3  (    '♦SCOi    360Ci  '♦or) 

YARDING   COST   ■   $  16392.63 

Lancing  building  rosT  -  $ 

HALLING  COST  -  S  AO.?? 
TCTaL   CCST   -  $  ?\7M7,B6 


lM  T   C  I£T4NCf 
*3        'tc76  FT 
26'47  FT 
i»l        139?  FT 
«?       aS'tS  FT 

F  .  CC 


KUh-BER   OF  LMTS  - 

kupber  of  landings  •  ? 
CASE  *  i: 


landing     coord INATFS ( X< Y J  7 ) 
#1  (    3900/    3100i  3001 


#2  (    5989/    7:35Mt  6971 

YARDING   COST   ■   $  15372.19 

Landing  building  cost  ■  s 

HALLING  COST  «  $  41.45 
TCTaL   cost   -  S  186P3.65 


LNIT  DISTANCE 


#1 
#2 
*3 


2368  FT 

2917  FT 

2598  FT 

2263  FT 


3P1C.CC 


CASE  « 


LANDING 
1i\ 


COORD InaTES ( X<  Y j  7 ) 
(    3900/    3100/  300) 


#2  (    4500/    3600/  400) 

YARDING  COST  ■  $     18919. 9C 

Lancing  building  cost  ■  s 

HALLING  cost  ■  •  31.85 
TCTaL    cost   -  S  21811.75 


LMT  CISTANCE 


#1 
#3 
#4 

#2 


2368  FT 

25S8  FT 

4259  FT 

2343  FT 


286C  .CO 


CASE   #  3: 


LANDING 
#1 


COORDINATFS ( V / Y>  7 1 
(    4437/    3548/  389) 


#2  (    1400/    2000/  200) 

yarding  cost  -  $  17504.07 
Landing  buIlding  cost  .  « 
halling  cost  ■  «  31.33 

TCTaL   CCST   -  e  20695.40 


LNIT  CISTANCE 


#2 
#3 
#4 
#1 


2397  FT 

294C  FT 

3854  FT 

1392  FT 


31f C.CC 


CASE   «  4! 


LANDING 
#1 


COORD INATFS ( X/ Yi 7  ) 
1    4437/    3548/  389) 


•2  (    3900/    3100/  300) 

yarding  cost  -  s  18813.47 
Landing  building  cost  -  « 
halling  cost  ■  s  33.96 

TCTaL   CCST   «  S  21647.43 


L^.  IT  CISTANCE 


#1 
*? 
#4 
#3 


3C46  FT 

2397  FT 

3854  FT 

2598  FT 


280C  .CO 


CASE  # 


LANDING     COORDINATFS ( X/ Yi 7 ) 
#1  (    4437/    3548/  389) 


#2 


(    5600/    5700/  600) 


LNIT  CISTANCE 

#1  3C46  FT 

«?  2397  FT 

#3  4276  FT 

#4  2647  FT 
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YARDINf,  COST   ■>  S  l(»2l3i>»f 
L4^CI^G   euILDlNR  COST   -  «  S^Sc.CC 
HAULING    COST    «    «  43.57 
TCTaL   CCST   -  S  217!i?.43 

CASE  «  6: 


LANDING 

CCCRO InaTF? ( X»  Y, 7  I 

LMT 

r;isTANCE 

n\  ( 

?0OC<  ?0C) 

♦  1 

139?  FT 

«3 

25CS  FT 

#2  ( 

5600'    5700*  600) 

#? 

2773  FT 

»>t 

2647  FT 

YaRC ING 

COST   ■   S  I'tSfiLai 

LANCING 

BUILDING  COST  ■  « 

CO 

HALL iNG 

COST  -  «           36. 8C 

TCTaL    CCST   ■  «  19018.11 


r*sE  D  T. 


LANtJlNC 

CCORD INATFS ( X<  Yi  7  1 

LN  IT 

C ISTANCE 

#1  ( 

4%37>    3548*  389) 

#1 

3046  FT 

#4 

3854  FT 

#2  ( 

4437*    3548*  389) 

»2 

2397  FT 

#3 

2940  FT 

YARR ING 

COST   ■   t  19240.34 

Lancing 

BUIL  DING  COST  ■  % 

28CC 

CO 

HALL ING 

COST    ■   $  34.80 

TCTaL   CCST   -  S  22075.14 


NUh-BFR  OF   LMTS   ■  4 
NUMBER     CF   LANDINGS  -  1 

CASE  #  i; 


lanoing 

C00RDlNATES(X<Y*7) 

LNIT 

CTSTANCE 

*\  ( 

4437'    3548*  389) 

#1 

3C46  FT 

*2 

2397  FT 

#3 

2940  FT 

*4 

3854  FT 

yarcing 

COST  •  «  19240.34 

Lancing 

BUILDING  COST  -  • 

140C 

.CO 

Hauling  cost  ■  «        34. 80 

TCTaL   CCST  ■  c  20675.14 


»  LANDING  location  ^'CCEL  ♦ 
»  BEST  SOLUTION  t 

♦  » 

««««««**¥*«¥¥¥¥¥¥¥lri*¥«*¥¥«*«««*«*«[l¥¥««¥¥¥««*««V¥1'«¥V****V¥* 

NO  CF  uMTS  "4  W 
NO  CF  LANDINGS  -     3  ¥ 

¥ 

»¥»¥¥¥¥*¥¥¥»¥¥¥»¥¥¥♦¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥»¥¥¥  ¥ 


LANDING  C00RDINATES(X/Yi7) 
#1  (    1400*    2000*  200) 

#2  (    4500*    3600*  400) 

#3  (    5989*    7354*  697) 

Yarding  COST  -  •  13327.89 
lancing  ruilding  cost  -  « 
Hauling  cost  -  •  36.61 

TCTaL  CCST   -  t  18584.50 


UNIT 

DISTANCE 

¥ 

«\ 

1392  FT 

¥ 

it  3 

2505  FT 

¥ 

«2 

2343  FT 

¥ 

»4 

2263  FT 

¥ 

5220  .CO 


¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥♦¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥»¥¥¥»¥¥»»»¥¥¥¥¥¥¥¥¥¥¥¥¥¥»¥¥ 
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HELICOPTER  LANDING  LOCATION  MODEL 
EXAMPLE   PROBLEM  #1 


Figure  8. — Plotter  output,  helicopter  logging  problem  1. 

18 


Problem  2 


In  problem  1,  the  variable  NALT  (columns  79-80,  card  2)  was  set  to  0.    Thus  the 
program  considered  locations  along  the  entire  road  in  100- foot  increments  (variable 
XINC,  column  67-76,  card  2)  for  potential  landing  locations.     In  problem  2,  let  us 
suppose  that  the  analyst  wished  to  consider  only  five  specific  locations  on  the  road 
as  possible  alternatives.    The  alternative  locations,  and  their  respective  landing 
building  costs  are  given  in  table  1.    Now  evaluate  several  different  allocation 
schemes;  first,  various  allocation  schemes  utilizing  two  landings,  and  also  four 
schemes  utilizing  three  landings.    These  schemes  are  shown  in  table  2.     Input  format 
for  this  problem  is  illustrated  on  a  FORTRAN  coding  sheet  in  figure  9.     It  should  be 
noted  that  NALT  (column  80,  card  2)  is  now  set  equal  to  5.    Also  note  that  landing 
costs  do  not  have  to  be  included  on  roadway  cards  since  these  costs  are  now  associated 
with  specific  landings.    Output  of  example  problem  2  is  given  in  figures  10  and  11  for 
the  line  printer  and  plotter,  respectively. 


Table  I . --Alternative  landing  data 


Alternative 

Coordinates 

Fixed  Cost 

Variable  Cost 

1 

(1,400;  2,000;  200) 

1,600 

0.0002 

2 

(4,500;  3,600;  400) 

1,100 

.0003 

3 

(5,600;  5,700;  600) 

1,300 

.0003 

4 

(6,000;  7,400;  700) 

1,300 

.0003 

5 

(6,450;  7,850;  750) 

1,300 

.0003 

Table  2. --Specific  allocation  schemes  with  three  landings 


Scheme  Landing  Unit 


1 

1,3 

3 

4 

5 

2 

2 

1,3 

4 

4 

5 

2 

1 

1 

2 

2 

3 

3,4 

1 

3 

2 

1,2 

4 

4 
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Figure  9. — Input  data^  heliaopter  logging  problem  2. 

20 


Figure  10. — Line  printer  output ^  helicopter  logging  problem  2. 
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Figure  11. — Plotter  output^  heliaopter  logging  problem  2. 
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Limitations 


Limitations  of  the  program  with  respect  to  the  number  of  critical  parameters  are: 

Vavametev  Maximum 

Number  of  units 
Number  of  road  segments 
Number  of  restricted  road  segments 
Ridges 


8 
20 
10 

1/unit 


The  first  three  limitations  shown  are  imposed  primarily  to  conserve  storage  require- 
ments and  can  be  easily  increased.     The  fourth  limitation  is  the  number  of  ridge 
constraints  that  can  be  considered  per  unit.     Presently,  the  program  is  limited  to, 
at  most,  one  ridge  constraint  per  unit  (units  may  have  a  common  constraining  ridge) . 
This  limitation  cannot  be  altered  easily,  as  a  change  would  entail  additional  program- 
ing.    Present  capabilities  will  prove  adequate  for  most  helicopter  logging  sales. 

When  evaluating  specific  alternative  landings  (NALT  >0) ,  hauling  costs  are  not 
calculated.     Under  this  option,  landings  need  not  be  located  on  an  existing  road, 
therefore  actual  hauling  costs  may  be  unknown  until  the  location  and  cost  are  speci- 
fied for  a  spur  road  from  the  haul  road  to  the  landing.     As  shown  by  the  preceding 
examples,  however,  these  costs  generally  prove  to  be  very  small  as  compared  to 
landing  and  yarding  costs . 
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Headquarters  for  the  Intermountain  Forest  and 
Range  Experiment  Station  are  in  C^den,  Utah. 
Field  Research  Work  Units  are  maintained  in: 

Boise,  Idaho 

Bozeman,    Montana  (in  cooperation  with 

Montana  State  University) 
Logan,    Utah  (in  cooperation  with  Utah 

State  University) 
Missoula,   Montana  (in  cooperation  with 

University  of  Montana) 
Moscow,  Idaho  (in  cooperation  with  the 

University  of  Idaho) 
Provo,  Utah  (in  cooperation  with  Brigham 

Young  University) 
Reno,   Nevada  (in  cooperation  with  the 

University  of  Nevada) 


